Rapid anomalous baryon number violating 'sphaleron' interactions present at high 
temperatures in the standard model certainly play an important role in determining 
the baryon asymmetry of the universe. Despite the large effort invested [1] however, 
we still lack a convincing theory of baryogenesis at the electroweak phase transition 
(EPT). What is much better established, is that baryon number was badly broken at 
temperatures above the EPT [2,3], which is known to greatly alter the more conventional 
pictures of baryon number production, e.g. in GUTs [4]. Moreover, the necessity 
of sphaleron reprocessing of particle asymmetries has made it possible to devise new 
schemes of baryogenesis [5-7]. For example in leptogenesis [5] one first generates an 
asymmetry in the lepton number L, which then becomes (partly) reprocessed into a 
baryon asymmetry by sphalerons. 

Loosely speaking the effect of sphaleron interactions is to break B + L while 
conserving B — L. Then the equilibrium conditions due to various 'ordinary' and 
anomalous interactions imply [3,8], to lowest order, that the B and L asymmetries are 
CN ' proportional to B — L. Thus a universe that initially has B — L = 0, as is the case 

in : 

in the simplest GUTs, should have zero baryon asymmetry. This realization added 
motivation for studying the generation of a baryon number at the EPT. Also, any source 
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of new exotic baryon or lepton number violation together with the sphalerons has the 
potential to destroy the baryon asymmetry, even with a nonzero primordial component 
of B — L. To avoid this outcome, it appeared to be necessary to place stringent limits 
on various effective operators [8-11] . These limits however, were shown to be weakened 

•i-H . 

^ ■ significantly due to the effective conservation of e^-number at temperatures T > 10 

TeV, and to only apply to flavour independent couplings [7,10,12,13] . 

In more accurate treatment of the equilibrium conditions [13,14], it was shown that 
even in the case B — L = - assuming that the EPT is weakly enough first order - 
the vacuum mass effects in the broken phase (re)generate a small but nonvanishing 
equilibrium value for B. Still it was maintained that above the EPT temperature Tc, 
B should be zero. In that case, should the EPT be of first order, as might yet be the 
case even in the minimal standard model [15], the baryon number would again be zero. 
(Unless of course, it was generated at the phase transition) 

In this letter we will show that the baryon and lepton numbers remain nonvanishing 
throughout the evolution of the universe, even in the case B — L = 0, if there are lepton 
flavour asymmetries. Above Tq a small nonzero B is preserved due to the nonidentical 
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dispersion relations of particles in a plasma. We will first find the effective thermal 
masses of particles in a primeval plasma, and then show how their presence alters the 
chemical equilibrium conditions so as to give rise to a nonzero B. We discuss the 
implications of our results to the Affleck-Dine [6] type baryogenesis, where large initial 
asymmetries may be generated [16]. 

When a particle propagates in a plasma, it continuously interacts with other 
particles in the background. Due to these interactions, for example the Dirac equation 
for a massless fermion becomes 



where a and b are some perturbatively calculable momentum and energy dependent 
functions [17] that vanish at the zero temperature limit, and u is the four velocity of 
the plasma. The operator multiplying \1/ in (1) is just the inverse propagator, the zeros 
of which define the dispersion relation E = E(k). The functional form of the dispersion 
relation is in general rather complicated, but one can introduce effective 'thermal' 
masses M(T) such that to very good accuracy [17] 



Thermal masses get contributions from both gauge and Higgs interactions and in general 
are different for different chiralities and families. Defining X{ = Mj(T)/T, one finds the 
following expressions in the standard model: 



where g,g' and g s are the usual SU(2), U(l)y and 577(3) couplings and the 
Yukawa couplings derive from the lagrangean Cy = J2ij(hij^i^ e jR + hfjQi&UjR + 
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[(l + a)^ + bi]^> = 0, 



(1) 



E(k) ~ ^k 2 + M 2 (T). 



(2) 




(3) 



hfjQi&djR) + h.c, with $ = ia^^- We choose the RH fermion basis such that h^h? is 
diagonal, the LH lepton basis such that h e h e ^ is diagonal, and the LH quarks such that 
h u h u ^ + h d h d ^ is diagonal. In this basis the eigenvalues of h^h and hh) are the same, 
and equal to 2m 2 /i> 2 , where v is the Higgs vacuum expectation value (= 246 GeV). 
Note that we do not know the eigenvalues of the matrix h u h u ^ + h d h d ^ . However, since 
the CKM mixing angles are small, we can assume that the eigenvalues ~ m 2 . + m 2 ,. . In 
any case, this is not important for our calculation. Numerically one finds 1 for example 
X£^l — 0.21, xg } R ~ 0.13 and x t L — 0.58 (see also table 1 below). 

Dispersion relations for bosons can be similarly derived and approximated by 
thermal masses. The thermal mass of the Higgs boson turns out to give xh — 0.59, 
assuming T 3> Tq, mn = 60GeV and m top = 174GeV; for a complete expression see e.g. 
ref. [12], but the actual value of Xh is not important here. The thermal masses of gauge 
bosons will turn out to be irrelevant for our purposes. 

Thermal masses can have an effect on baryon number through the boundary 
conditions, such as charge neutrality Q em = (cf. equation (7) below). Total charges 
are proportional to asymmetries in particle densities, whereas the equilibrium conditions 
are relations between chemical potentials. It is then clear that finite mass effects can 
modify the solutions to equilibrium equations subject to boundary conditions such that 
nontrivial solutions may exist [14]. 

Indeed, the occupation number of a given particle species with momentum k (in 
the rest frame of the heat bath) is given by 

n(k, T, fn) = g(eM(E(k) - m)/T) ± (4) 

where +(-) refers to fermions (bosons), and g is the number of internal degrees of 
freedom of the particle. Then the asymmetry in that particle species is 

/d 3 A; 
— -(n(/c,T, m) -n(k,T,-Vi)) 
( 2y (5) 

= ^T 2 )x(a±-^), 

where we have expanded to first order in ^, and a + = 1 and a_ = 2. Assuming that 
the dispersion relation is of the form (2), the function 5 becomes (In the notation of 



1 We assumed that the coupling constants take their values measured at the weak scale. More 
accurately we could write, e.g. x\ L ~ 0.23% + 0.10, with \ = ln(Mz / Aqcd) / ^ n (T / Aqcd) ■ 



ref. [13] 5 = a± — a.) 
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6± = a± --zJ x dyy ^^J^TW (6) 

For fermions one finds to first order ~ ^f^^ 2 ' wmcn * s a ver y good approximation for 
all quarks and leptons. In the bosonic case the integral expression for 5- is nonanalytic 
near x = 0, but one can show that for x < 1 to good accuracy 5- ~ x. 
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Table 1. Shown are the effective masses x = M(T)/T and corresponding values of 5-functions 
as given by equation (6). We took the QCD coupling corresponding to the scale Mz (see footnote 
1). Only the top quark Yukawa coupling is large enough to show up in these results. We have 
used mtop = l74GeV and a s = 0.11. 

We will now solve for the baryon number in the presence of finite temperature 
mass effects in a universe with B — L = 0, but Li ^ 0. At T > Tc, there are 
additional boundary conditions of vanishing total charge- and isospin densities, so that 
the complete set of boundary conditions to be used is: 

Qem = Qs = B-L = 0. (7) 

Rapid decays and inverse decays (when kinematically allowed) and various 2 — 2 
scattering processes induce the usual set 2 of 'ordinary' equilibrium relations between 

2 We can assume that the temperature is less than (D(few)TeV, so that also right 
handed electrons are in thermal equilibrium [12]. 
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chemical potentials [8]: 



(1- + HO = fJ>u R — Hu L = ^0 

^d R - Vd L = ~llQ fJ>i R - Vi L = -V0 ( 8 ) 

Hd L — Hu L = V-W V>i L — fl>i = /J>W, 

where the subscipts indicate particle species - - "0" and " — " being the Higgs doublet, 
% a neutrino species, and il and i R are charged leptons of the ith generation. 
Additionally, there is an equilibrium condition due to anomalous sphaleron processes 
which, with help of (8), can be written as 

9fJ, UL + 6/j w + A* = 0, (9) 

where fx = £V We can use the equations (5) for particle asymmetries along with the 
equilibrium conditions (8) to write the global charges as 

Q = (6 - 2A U + A d ) tiuL - (18 - 45 w - 25- - A d - A e ) m 



(10) 



+ (14 - 2A uR - A dR - A eR - 25-) /j -2/j + A/j l + A^ R 

3 1 

Q 3 = -(11 - A5 W -5- - -A dL - - A eL ) fj, w 

3 1 

+ -^(Acil - A uL ) [i uL + -(A^ L - A//„) + fj, (5 - 5-) 

B = (12 — A q ) h u l + (6 — A d ) [i w + (A dR — A uR ) fx Q 

L = S/j, - Ay, v - A/j l - A/j r + (6 - A e ) fi w - (3 - A eR ) fj, , 

where we have dropped the common multiplicative factor A Jr g T , (we define L = An/s, 
where s is the entropy density so that g* is the number of relativistic degrees of freedom 
taken to be 106.75 from now on) and defined 



Afi v 


= ^2/Ji5i 

i 




A/jx 


i 


X = L,R 


A/x 


= 6 f^ 


X = L,R fi = Ui,di,ei 



(11) 



and moreover Af = AfL + Af R and A q = A u + A d . Notice that the chemical potential 
Hi in definitions (11) is always that of neutrinos. We must now apply the boundary 
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conditions (7) to the expressions (10). At first sight this looks rather messy, but 
fortunately some of the coefficients in (10) turn out to be zero, so that an analytic 
solution of decent length can be obtained for B. 

First of all, dispersion relations of particles within the same weak doublet are 
identical, so that, Sf uL — Sf dL = 0, So — 5- = 0, and AfiL = Afi v . Then only the term 
proportional to fiw remains in the expression for Q% and therefore, even with thermal 
corrections we get fiw = 0. After this simplification it is rather straightforward to 
proceed: we use (9) to solve fi u L = — /u/9 and then solve fio from the charge equation: 

fM) = u ^ Ai {(| - A 2^ " A ^L ~ A// fl j , (12) 

where we used the shorthand notations Ai = 2A uR + A dR + A eL + 2S and 
A 2 = |(2A U — Ad). Finally, the last constraint B — L = becomes a consistency 
condition for fx: 

13 A ? / 8 A ^ 3 + A3 1 , ^ 3 + A 3 



3 9 - ( 3" A2) 14^) ^ = A^ + (1- II -^)(A^ + A^), (13) 

where A 3 = A^r — A u r — A e R. If all thermal corrections vanish, equation (13) has only 
the trivial solution fj, = 0, which would immediately lead to B = L = 0, as indeed was 
earlier thought to be the case [8,13,14] . 

At first sight, (13) might be thought to be giving rise to relatively large 
asymmetries, because the gauge contributions to effective masses are not particularily 
small. However, the gauge contributions are the same for all families and their added 
contribution is proportional to fx. Indeed, we can write the A//x-terms in the R.H.S of 
equation (13) as 

&H X = 8 9 x n + Y,8r x n i , X = L,R (14) 

i 

where S 9 and S Y respectively are the gauge and Yukawa mass corrections (see equation 
(3)). Thus the equation (13) is of the form A\i = afi + c, where c is nonzero, and hence 
a nontrivial solution exists, only because lepton family number is conserved and the 
lepton Yukawa interactions differ from one family to another. This is in close analogy 
to the situation when T <Tc [13,14]. Using (14) one can rewrite (13) as 



^ ( S A -, 3 + A 3 3 + A 3 g g \ 

(y--Q-- ( 3- A2) 143A;-^-( 1 -i43A^)^ + ^)j^ 



(15) 



which shows explicitly the dependence on Yukawa terms A\i^ = ^i^A*i- I n principle 
it is easy to trace backwards the steps from (15) to (12) to (10) and solve for B exactly. 
However, glancing at the size of the various thermal correction terms (see also table 1), 

A q ~ 0.49 Ai ~ 1.39 A 2 ~ 0.03 A 3 ~ -0.03 5 9 L ~ 0.007 6 9 R ~ 0.002, (16) 



one can see that all of the quantities in (16) are small enough to be neglected. Then 



one finds a simple expression for B oc 12fi UL ~ — ifj,, 



B.-^AkM). (17) 

Including all A-terms (with Xh = 0.59) one obtains B ~ —4.15 x 10~ 3 A / u^/T, which is 
less than 3 percent off from the value given by (17)! The Yukawa couplings of leptons 
are rather small and obey the hierarchy h e <C <C h T (see table (2)), so that to very 
high accuracy A(j?[ ~ 5 Y /j t oc S^L T . 
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Table 2. Shown are the Yukawa contributions to the lepton thermal masses and the correspond- 
ing values of 5 functions. For right handed charged leptons use x\ = 2x\ and 5r = 25 l- 

We still need to express B in terms of primordial quantities. It is easy to show 
that if B — L = and |S — Lj are conserved, then 

L T (t) = ^([B + L](t) + 2AL Te + 2AL TfI ). (18) 

where L T = L TL + L TR + L Ut and ALij = L; — Lj. The sphalerons force B + L and 
the Higgs chemical potential to be zero, up to the small mass corrections that we are 
calculating, so that L T ~ A ^ g and 

B * -^^ ^q^ - L3 >< 10 " 7 ( AL - + A V) , (19) 
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This is our final result. It should be noted that mass effects alone are not enough 
to provide a nonzero B, although they set the scale of B in terms of primordial 
asymmetries. Instead, primordial values of leptonic asymmetries must be different, i.e. 
ALij ^ 0! 

Let us compare our result (19) to that of refs. [13,14], where the situation at 
T < Tq is considered. Equation (23) for B in ref. [13] is similar to our equation (17) 
and it can also be brought into the form (19). Neglecting all mass corrections except 
the bare lepton Yukawa terms responsible for a nonzero result, one obtains from ref. 
[13] B ~ few x 10~ 7 (AL er + AL MT ). Of course, if the EPT is of second order or weakly 
enough first order, it is the latter result that is relevant. However, if the EPT is of first 
order, then sphaleron processes drop out of equilibrium instantaneously in the broken 
phase, the equilibrium conditions below Tq never get realized and the value of B given 
by (19) gets frozen into the system. The point we wish to make is that either way, a 
nonzero fraction of primordial baryon number survives. 

In order for a component of primordial baryon number to survive at the level of 
the observed asymmetry B now = ns/s ^4 x 10 -11 [18], primordial asymmetries must 
be of order ~ 10 -4 . This is perhaps at best marginally possible in the simplest GUT 
scenarios, but it is much more feasible in the Affleck-Dine mechanism [6], where a 
baryon asymmetry is generated during the decay of a scalar condensate, which can 
have a large component of baryon or lepton number stored within. It was previously 
found [16] that sufficiently large (B > 10~ 2 ) primordial B + L asymmetries can survive 
in the Affleck-Dine scenario, because the squark condensate does not evaporate before 
the electroweak phase transition, giving the W a large mass and keeping the sphalerons 
out of thermal equilibrium. However, a large amount of entropy production is then 
required after the phase transition. We have shown here that within the Standard 
Model, given a slightly smaller "large" primordial asymmetry ~ 10~ 4 , it would not be 
washed out, but rather be diluted by the sphalerons to give approximately the right 
baryon asymmetry today, so that no entropy production is required. 

Let us finally discuss the apparent limitation that nonzero ALj/s are required for 
the mechanism to work. First of all, it is not at all unnatural to obtain differences in 
the asymmetries that are of the order of the asymmetries themselves. On the contrary, 
in e.g. A-D models it is natural to expect that the scalar condensate initially lies in 
some arbitrary direction in the slepton space, which in general is not symmetric in 
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flavours. In such case the subsequent evolution and decay of the condensate naturally 
gives rise to large differences in lepton asymmetries. Thus large ALy's certainly can 
be produced. Secondly, in models with B — L symmetry ALy s are always conserved 
quantities. Moreover, even in the case that the lepton number was broken at some 
higher temperature scale due to some new exotic interactions [9], B typically would 
not vanish. This is due to the approximate e^-number conservation, which protects 
the right handed electron asymmetry, and hence a nonzero AL Te at any temperature 
T > 10 TeV [12]. Thus, except in some very speculative cases (for details see ref. [12]) 
the survival of a nonzero B is always guaranteed. 

In conclusion, we have computed the equilibrium value of the baryon asymmetry B 
in the early universe above the electroweak phase transition temperature Tq. We have 
found that due to thermal effects (nontrivial dispersion relations) a small, but nonzero 
value of B persists throughout the evolution of the universe, even in the case B — L = 0. 
The final value of B is roughly seven orders of magnitudes smaller than the primordial 
value, with the scale set by the Yukawa coupling of the tau lepton. Our results 
complement the analysis of refs. [13,14], where the vacuum mass effects at the broken 
phase at T < Tc were shown to also secure a nonzero B. Combining these results one 
has the proof of persistence of a fraction of primordial baryon asymmetry regardless of 
the order of the electroweak phase transition. This result favours Affleck-Dine type 
baryogenesis models, which otherwise could produce too large a final baryon asymmetry; 
the scale of relaxation of B due to modified equilibrium processes can make these 
theories work even for rather large initial asymmetries. 
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Abstract 

We consider the evolution of baryon number B in the early universe under the influence 
of rapid sphaleron interactions and show that B will remain nonzero at all times even 
in the case of B — L = 0. This result arises due to thermal Yukawa interactions that 
cause nonidentical dispersion relations (thermal masses) for different lepton families. 
We point out the relevance of our result to the Affleck-Dine type baryogenesis. 



